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ABSTRACT 


An  experimental  shock  tube  program  was  undertaken  in  order  to  calibrate  electro¬ 
static  probes  which  may  be  used  to  make  point  ion  density  measurements  in  continuum 
flowing  plasmas. 

The  empirically  derived  formula  which  relates  the  ion  density  to  the  measured  probe 
current  was  obtained  in  terms  of  the  appropriate  similarity  parameters  and,  over 
the  range  of  conditions  covered,  was  found  (for  spherical  probes)  to  be  expressible 
in  the  form: 

1 

1  +0.  61  t,n  (1+r/A-J 

i  =  (-0)  U  [  1.  2  *  0. 09  Re] 

where  i  is  the  Nusselt  number  for  mass  transfer  of  the  positive  ions,  (-0)  is  a 
dimensionless  probe  potential,  X^/r  is  the  ratio  of  the  Debye  length  to  probe  radius, 
and  Re  is  the  flow  Reynolds  number  based  on  the  probe  diameter.  This  empirical 
correlation  was  used  to  interpret  the  ion  currents  collected  by  a  radial  array  of 
negatively  biased  spherical  probes  which  were  positioned  so  as  to  Intercept  the  ionized 
wake  of  a  hypervelocity  model  in  a  ballistics  range.  The  integrated  ion  densities  were 
found  to  be  in  good  accord  with  the  integrated  electron  densities  which  were  simul¬ 
taneously  measured  using  microwave  methods.  The  application  of  the  electrostatic 
probe  measuring  technique  to  the  study  of  negative-ion  chemistry  in  wakes  is  also 
illustrated  and  briefly  discussed. 
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I.  INTRODUCTION 


The  desirability  of  being  able  to  make  point  electron  density  measurements  in  the 
turbulent  wakes  of  hypervelocity  projectiles  has  long  been  recognized  as  a  necessary 
step  toward  gaining  a  proper  understanding  of  their  structure  and  their  radar  scatter¬ 
ing  characteristics.  Although  ballistics  ranges  do  provide  an  appropriate  laboratory 

environment  for  making  such  measurements  and,  Indeed,  measurements  of  average 

(1  21 

electron  densities  as  a  function  of  distance  behind  the  projectile  have  been  obtained,  ’ 
no  reliable  point  measurements  of  electron  density  fluctuations  have  so  far  been  made. 

In  a  recent  paper,  Sutton'  'proposed  a  ballistics  range  experiment  which  might  allow 
such  measurements  to  be  performed  with  cylindrical  collision-free  Langmuir  probes, 
held  in  tension  so  as  to  avoid  the  problem  of  probe  bending.  In  order  to  operate  the 
probe  in  the  collision-free  regime  and  at  the  same  time  have  a  turbulent  wake, 
relatively  large  ballistic  models  and  low  environmental  pressures  are  required. 
(Sutton's  calculations  were  based  on  the  theoretically  computed  wake  conditions  be¬ 
hind  a  6.  8  cm-diameter  sphere  traveling  at  4.  72  km/sec  in  an  ambient  pressure  of 
10  torr. )  Unfortunately,  the  ballistics  range  experiments  carried  out  to  date  at 
AC  Electronics  -  Defense  Research  Laboratories'  (AC-DRL)  facility  have  been  con¬ 
ducted  with  smaller  models  (typically  spheres  which  range  in  diameter  from  0.  5  cm 
to  2.  5  cm)  and  at  pressures  greater  than  10  torr. 

There  are  a  variety  of  reasons  which  lead  to  the  choice  of  a  particular  set  of  condi¬ 
tions  for  operation  of  a  ballistics  range.  These  depend  upon  the  type  of  measurement 
that  is  to  be  made;  however,  one  of  the  most  Important  factors  which  predicates 
against  large-scale  models  is  the  fairly  obvious  economic  one  that  the  costs  of  the 
ballistics  range  facility  and  its  operation  Increase  at  a  greater  than  linear  rate  with 
increasing  model  size  for  a  typical  reentry  speed  of  about  6  km/sec.  It  seemed 
reasonable  therefore,  to  examine  the  possibility  of  using  electrostatic  probes  which 
can  be  operated  in  the  continuum  regime,  where  the  mean-free-path  is  smaller  than 
the  characteristic  probe  size.  The  obvious  advantage  to  be  gained  by  adopting  this 
approach  is  that  one  can  continue  to  use  the  existing  ballistics  ranges  and  their 
established  methods  of  operation  to  investigate  wake  structure  ever  a  considerable 
range  of  body  sizes  (over  at  least  a  factor  of  5  in  scale  size). 
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Clearly,  the  crucial  problem  in  using  continuum  electrostatic  probes  in  the  range 
environment  is  the  method  of  interpreting  the  probe  current  measurements.  Although 
a  number  of  theoretical  treatments  of  probe  response  in  continuum,  flowing  plasmas 
are  available  (e.  g. ,  Lam,^  Hoult,^  Sonin,1 ^  and  Su.^to  mention  only  a  few  repre¬ 
sentative  works),  the  complexity  of  the  situation  is  such  that  each  analysis  is  addressed 
to  specific  limiting  conditions  of  flow,  probe  geometry,  etc. ,  and  a  unified  theory 
spanning  the  wide  range  of  conditions  encountered  in  a  ballistics  range  is  not  at  present 

available.  Consider,  for  example,  a  probe  about  one  millimeter  in  diameter  under 

(2) 

typical  range  conditions.  Except  in  the  very  near  wake  (less  than  about  a  hundred 

(81 

body  diameters  downstream  of  the  projectile)  the  flow  in  the  wake  is  subsonic,  and 

the  primary  parameters  governing  the  probe  response  are  expected  to  be  the  local 

Reynolds  number  based  on  probe  diameter,  Re,  the  ratio  of  Debye  length  to  probe 

radius,  X^/r,  and  the  dimensionless  probe  potential,  0  =  eV/kT.  As  the  probe 

passes  from  the  near  wake  to  the  far  wake,  the  Reynolds  number  ranges  from  about 

-2 

200  to  4,  and  the  Debye  length  ratio  from  about  10  to  1,  as  illustrated  in  Figure  1. 

The  available  theories,  on  the  other  hand,  are  tailored  largely  for  the  limiting  condi¬ 
tions  of  either  very  small  or  very  large  Debye  lengths  and  very  small  or  large 
Reynolds  numbers,  as  indicated  in  the  figure  (note  that  the  situation  in  Figure  1  is 
somewhat  oversimplified,  since  the  third  coordinate  of  the  similarity  parameter  space, 
the  probe  potential  eV/kT,  has  been  omitted).  Another  deficiency  in  the  available 
theories  is  their  inability  to  deal  with  ion  collection  by  the  downstream  side  of  the  probe, 
where  flow  separation  occurs;  thus  in  a  turbulent  wake  it  is  generally  impracticable 
to  use  probes  with  collecting  surfaces  always  oriented  into  the  flow. 

It  is  apparent,  therefore,  that  an  experimental  calibration  of  probes  is  required  in 

flow  conditions  simulating  those  of  the  ballistics  range.  Such  a  calibration,  and  its 

application  in  an  investigation  of  hypersonic  wakes,  is  presented  in  this  report.  The 

probe  calibration  experiments  were  carried  out  using  a  shock  tube  where  the  variation 

of  the  electron  density  behind  the  primary  shock  could  be  measured  using  a  focused 

(2  9) 

microwave  interferometer.  ’  '  Most  experiments  reported  here  were  done  with 
spherical  probes  so  as  to  eliminate  the  Influence  of  flow  direction  on  probe  response, 
and  only  ion  collection  at  negative  probe  potentials  wa>»  considered.  The  range  of 
the  probe  Reynolds  number  and  the  ratio  of  the  Debye  length  to  probe  radius  covered 
in  this  calibration  is  also  illustrated  in  Figure  1. 


Figure  1  Range  of  Probe  Re  andXD/r  Covered  in  a  Typical  Wake  and 
also  in  the  Present  Calibration.  Also  shown  are  limits  of 
available  theories. 


The  spherical  probe  calibration  part  of  the  experiments  described  here  is  discussed 
in  Section  II  (see  Appendix  A  for  cylindrical  probe  calibration)  and  its  application  to 
measuring  mean  ion  density  profiles  in  the  wakes  of  hypervelocity  projectiles  is 
demonstrated  in  Section  III.  Although  such  probes  are  also  capable  of  providing 
quantitative  information  on  the  fluctuating  components  of  the  ion  density,  this  report 
is  limited  to  a  discussion  of  the  mean  density  profiles. 
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II.  THE  SHOCK  TUBE  CALIBRATION 


A.  EXPERIMENTAL  APPARATUS 

A  pressure-driven  stainless -steel  shock  tube  having  an  internal  diameter  of  15  cm 
was  chosen  as  a  means  of  providing  a  flowing  plasma  whose  temperature,  pressure, 
and  flow  velocity  in  the  region  behind  the  primary  shock  could  be  accurately  calcu¬ 
lated  from  the  known  initial  conditions  and  measured  shock  velocity.  Details  of 
construction  and  the  method  of  operation  of  the  shock  tube  have  been  described  else¬ 
where  in  the  literature^ ^^and  will  not  be  repeated  here. 

The  test  gas  used  was  argon,  to  which  small  amounts  of  acetylene  and  oxygen  (1% 

C  ,H.,  and  5%  0>2)  were  added  to  produce  measureable  ionization  levels  behind  the 
shocks  at  relatively  low  temperatures  by  the  process  of  che mi- ionization.  The  pres¬ 
sure  behind  the  shock  (p2)  was  varied  between  5  and  1200  torr.  In  these  calibration 
experiments  it  was  not  necessary  that  the  precise  reaction  kinetics  of  the  chemi- 
ionization  process  should  be  known,  since  the  temporal  variation  of  the  electron 
density  was  measured  directly,  using  a  microwave  interferometer  (see  below). 

However,  it  is  generally  believed  that  the  dominant  positive  ion  which  persists  after 
the  initial  reaction  zone  in  lean  acetylene  flame^1*’  *^and  shock-tube  studies^1®’  ^ 
is  HgO+,  which  is  then  removed  via  a  dissociative  recombination  reaction  of  the  type: 

HgO+  +  e"  -  HgO  +  H  (or  other  neutral  products)  (1) 

Since  it  was  necessary  to  have  some  reasonably  reliable  estimate  of  the  diffusion 
coefficients  of  the  charged  species  for  the  method  of  data  reduction  used,  HgO+  was 
assumed  to  be  the  dominant  positive  ion  in  the  plasma  decay  region  behind  the  shock 
wave.  (Even  if  other  ions  such  as  CHO+  made  significant  contributions  to  the  total 
positive  ion  density,  it  is  believed  that  this  would  not  constitute  a  significant  error 
since  U  would  enter  only  through  the  diffusion  coefficient,  which  is  inversely  pro¬ 
portional  to  the  square  root  of  the  mass. )  The  only  negatively  charged  species  present 
were  electrons,  since  the  temperature  range  of  the  present  experiments  (2000°- 7000°K) 
was  too  high  to  allow  any  significant  concentration  of  molecular  negative  ions  such  as 
O'  and  02  to  be  formed. 
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The  experimental  arrangement  in  the  teat  section  of  the  shock  tube  is  illustrated  in 
Figure  2.  The  passage  of  the  primary  shock  wave  through  the  test  gas  produced  a 
region  of  extremely  rapid  ionization  which  was  followed  by  a  region  of  much  slower 
decay.  Within  these  regions  the  electron  density  was  measured  with  a  focused  micro- 
wave  probe,  ’  ’  and  the  ion  density  was  monitored  with  spherical  ion  probes. 

The  70GHz  focused -beam  microwave  interferometer  used  to  make  the  electron 
density  measurements  in  the  shock  tube  operates  on  essentially  the  same  principles 
as  reported  previously/®'  ^  The  main  virtues  of  this  system  are  that,  througn  the 
use  of  focusing,  high  spatial  and  temporal  resolution  is  achieved  In  the  plasma- 
interaction  region.  Also,  through  appropriate  design  there  is  essentially  free-space 
propagation  of  plane  waves  in  the  focal  region,1 ^permitting  very  straightforward 
interpretation  of  the  measured  quantities. 

A  schematic  view  of  the  microwave  module  is  shown  as  part  of  Figure  2.  Dielectric 
lenses  are  used  to  focus  the  microwave  energy  into  the  flow-through  section  which 
samples  the  near-axis  flow  in  the  shock  tube.  An  f/2  lens  was  chosen  to  give  a  depth 
of  field  greater  than  the  width  of  the  flow-through  section  (6. 4  cm  vs  3. 8  cm). 


Figure  2  Electrostatic  Probes  and  Focused  Microwave  Interferometer 
in  Shock  Tube 
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Measurements  with  Irises  lined  with  microwave  absorber  showed  that  90%  of  the 
energy  was  confined  In  a  cylinder  1. 6  cm  In  diameter.  Flush  dielectric  windows  are 
Inset  In  the  sides  of  this  section  to  permit  transmission  of  the  focused  beam  through 
the  ionized  flow.  The  thickness  of  the  windows  was  machined  to  an  Integral  number 
of  half-wavelengths  to  eliminate  multiple  reflections  of  the  microwave  energy.  The 
measured  standing-wave  ratio  In  the  box  was  less  than  1. 005.  The  diameter  of  the 
windows  (3. 8  cm)  was  sufficiently  larger  than  the  focused-beam  diameter  so  that 
there  were  negligible  effects  on  the  transmission  properties  of  the  focused  beam. 

The  receiver  circuitry  used  autodyne  mixing  with  enough  bias  for  the  mixer  provided 

by  the  reference  path  that  the  output  Is  A(t)slne(t),  where  A(t)  and  0(t)  are  the 

time-varying  amplitude  and  the  phase  change  of  the  signal  transmitted  through  the 

ionized  flow.  Through  the  use  of  an  additional  receiver  circuit,  A(t)  cos  6(t)  is 

also  recorded,  permitting  A(t)  and  9(t)  to  be  derived  explicitly.  These  quantities 

are  then  converted  to  collision  frequency  and  electron  density.  The  dynamic  range 

at  this  Instrument  Is  about  10^  to  10*3  cm"^;  the  electron  densities  In  the  tests 

were  in  the  range  of  1010  to  SxlO^cm'3  and  the  plasma  was  therefore  always  very 

underdense.  The  largest  value  of  electron  collision  frequency  measured  (for  p,*.  1000 

1 1  •  1  * 

torr  and  Re  «  2000  to  8000)  was  2x  10  sec  .  This  Is  still  well  below  the  operating 

frequency  (4. 4x10* 1  radians  sec"1).  Since  the  measured  attenuation  was  always  below 

(171 

5%,  the  phase  change  Is  almost  directly  proportional  to  the  electron  density.  ' 

Figure  3a  shows  a  typical  oscilloscope  display  of  the  A(t)  sin  0(t)  signal  measured 
behind  an  Incident  shock. 

The  spherical  ion  probes  (diameters  0. 76  mm  and  3. 2  mm)  were  mounted  on  an  In¬ 
sulated  sting  (0. 5  mm  In  diameter)  which  protruded  about  one  probe  diameter  in 
front  of  the  leading  edge  of  a  thin  wedge-shaped  support.  This  probe  assembly  was 
positioned  a  short  distance  ahead  at  the  microwave  box  and  sufficiently  far  off  the 
shock-tube  axle  that  no  disturbance  from  it  could  Influence  the  flow  within  the  box. 
Nearly  all  the  measurements  to  be  reported  were  made  using  this  configuration. 

The  probes  were  biased  at  voltages  ranging  from  -1. 4  to  -30  V  relative  to  the  shock- 
tube  ground.  A  block  diagram  of  the  electrostatic  probe  circuit  Is  shown  In  Figure  4. 
The  voltage  signal  generated  by  the  passage  of  the  probe  current  through  the  load 
resistance  was  passed  through  an  X100  operational  amplifier  and  the  resulting  voltage 
signal  was  displayed  on  an  oscilloscope.  This  system  thus  gave  both  large  bandwidth 
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Figure  3  Response  of  70GHz  Microwave  Interferometer  A  sinS(t)  and 
Electrostatic  Probes  to  Shock-Ionized  Acetelyne-Oxygenand 
Argon  Mixture 


and  dynamic  range.  Figure  3b  shows  the  measured  probe  currents,  displayed  on  a 
logarithmic  vertical  scale,  from  two  probes  biased  at  -15  V  and  -30  V;  the  displays 
shown  in  Figures  3a  and  3b  were  obtained  behind  an  Incident  shock  in  the  same  experi¬ 
mental  test.  One  major  division  on  the  vertical  scale  of  Figure  3b  corresponds  to 
approximately  a  factor  of  10  In  probe  current. 

A  few  measurements  were  made  behind  reflected  shocks  where  the  gas  velocity  is 
zero  (or  very  close  to  it).  For  these  measurements,  a  plug  was  placed  in  the  end  of 
the  centrally  located  flow-through  section  3.  8  cm  beyond  axis  of  the  focused  micro- 
wave  beam.  The  electrostatic  probes  were  mounted  on  the  plug  face  and  stood  1. 0  cm 
off  the  end  wall  (2. 8  cm  beyond  the  focused  microwave  beam  in  the  box).  For  data 
analysis  purposes  it  would  have  been  desirable  to  locate  both  probes  at  the  same  axial 
position  along  the  box  In  order  to  sample  elements  of  test  gas  with  identical  histories 
behind  the  incident  and  reflected  shocks  at  the  points  of  measurement.  This  was  not 
possible,  however,  because  the  electrostatic  probes  and  their  mounts  would  have 
seriously  interfered  with  the  focused  probe  measurements. 


Figure  4  Block  Diagram  of  Electrostatic  Probe  Electronics 
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Figure  5a  shows  oscilloscope  traces  (or  A(t)  sin  6(t)  and  A  cos  0(t)  (i.  e. ,  electron 
density)  taken  in  the  reflected  shock  configuration.  Figure  5b  shows  the  electro¬ 
static  probe  currents  (on  a  linear  vertical  scale)  from  two  probes  for  the  same  shot. 
The  lower  trace  is  for  a  spherical  probe  0.  76  mm  in  diameter  biased  at  -15  volts; 
the  upper  trace  is  for  a  cylindrical  probe.  The  undisturbed  gas  ahead  of  the  shock 
is  labeled  region  (1).  Region  (2)  is  the  gas  behind  the  incident  shock  which  shows 
the  rapid  chemi- ionization  and  the  slower  recombination  (Eq.  (1) ).  Region  (5)  is 
the  reflected  shock  region.  The  flow  velocity  of  the  gas  is  brought  to  zero;  the 
temperature  suddenly  rises  from  5, 150°K  (i.  e. ,  Tj)  to  9, 100°K  (1.  e. ,  Tg)  and 
chemi -ionization  starts  again,  sending  all  the  traces  off  scale.  The  electrostatic 
probes  follow  the  electron  density  quite  well. 

B.  SIMILARITY  PARAMETERS  AND  CORRELATION 

OF  SPHERICAL  PROBE  DATA 

In  order  to  calibrate  a  probe  in  conditions  which  do  not  precisely  duplicate  those 
where  it  is  to  be  applied,  it  is  necessary  to  establish  the  similarity  parameters 
which  govern  the  probe  response.  These  can  be  found  formally  by  dimensionless 
analysis.  To  facilitate  this,  the  following  observations  can  be  made  about  the  flow 
conditions  in  the  present  context: 

(1)  The  flow  about  the  probes  was  in  or  near  the  continuum  regime  in  the 
ballistics  range  as  well  as  in  the  shock  tube  where  the  calibration  was 
performed. 

(2)  Although  in  the  hypersonic  wake  the  flow  was  turbulent  (unsteady),  the 
flow  over  the  probes  was  assumed  to  be  quasi -steady,  which  would  be 
true  if  the  probe  diameter  were  small  compared  with  a  characteristic 
turbulent  eddy  size.  Steady  flow  was  established  in  the  shock  tube. 

(3)  The  flow  over  the  probes  was,  to  a  good  approximation,  chemically 
frozen  in  both  test  environments.  The  rate  coefficient  (a)  for  the 
dissociative  recombination  reaction  (1)  has  been  measured^)  over 
the  temperature  range  2,500°-5,000°K  to  be  a  -  1.  OST'^cm^sec"1. 

The  characteristic  plasma  decay  time,  defined  as  [d(ln  n ij/dt]"1*  1/an  , 
was  thus  found  to  be  somewhat  in  excess  of  100  ^sec  in  the  shock- 
tube  experiments.  Since  this  time  was  so  much  longer  than  the 
characteristic  flow  time  over  the  probes  (^  lpsec),  the  assumption 
of  a  chemically  frozen  flow  was  almost  certainly  Justified  even  after 
allowance  is  made  for  the  presence  of  a  weak  bowshock  wave. 

(4)  The  gas  was  weakly  ionized. 
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(b)  Electroatatic  Probes.  Top  trace  cylindrical  probe  (X  to  flow 
direction),  lower  trace  spherical  probe.  Note  different  time 
bases  between  9(a)  and  8(b) 


Figure  5  Response  of  Probes  in  Reflected  Shock  Region 
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Under  such  flow  conditions  the  current  drawn  by  a  probe  is  governed  by  seven 
equations:  the  (frozen)  ion  and  electron  conservation  equations,  Poisson's  equation, 
the  conservation  equations  for  mass,  momentum,  and  energy,  and  the  equations  of 
state  for  the  gas.  The  solution  of  these  equations  depends  on  the  following  physical 
quantities: 

r  =  probe  radius 

V  =  potential  of  probe  with  respect  to  plasma 

Vj  =  flow  velocity  of  gas 

f>  =  gas  density  (local  free-streain  values) 

u  =  gas  kinematic  viscosity 

T  =  gas  temperature 

nt,ne  =  ion  and  electron  number  densities  in  gas  (assumed  equal) 

D.,D  =  ion  and  electron  diffusion  coefficients 

l  e 

y  =  ratio  of  specific  heats  of  gas,  Cp/Cv 
Tw  =  probe  wall  temperature 

m  =  molecular  weight  of  gas 

k,cQ  =  Boltzmann's  constant,  premittivity  of  free  space 
e  =  electronic  charge. 


From  dimensional  analysis  we  find  that  eight  Independent  dimensionless  parameters 
can  be  formed  from  the  above  quantities: 

Re  =  Vj2r/u,  Reynolds  number  based  on  probe  diameter 

Xp/r  =  ratio  of  Debye  length  to  probe  radius,  where  Xp  -  ^e^n  /«  kT 

0  =  eV/kT,  dimensionless  probe  potential  with  respect  to  the 

ionized  gas 

M  =  vy^yRT  ,  Mach  number,  where  R  is  gas  constant  per  unit  mass 
y  -  ratio  of  gas  specific  heats 

Sc  «  u/Dj,  Schmidt  number  based  on  ion  diffusion  coefficient 

D./t)  =  ratio  of  ion  and  electron  diffusion  coefficients 

t  e 

Tw/T  =  dimensionless  probe  wall  temperature. 


It  follows  that  if  we  define  a  dimensionless  probe  current  i , 

i  Ir 

enjDj(4irr  ) 
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where 

I  =  total  current  drawn  by  probe, 

then  this  dimensionless  current  must  be  a  function  only  of  the  eight  dimensionless 
parameters  listed  above,  i.  e. , 

i  =  i(Re,  XD/r,  0,  M ,  y,  Sc  ,  ,  T^/T)  (3) 

We  note  that  i  is  simply  the  Nusselt  number  (Nu)  for  mass  transfer  (often  called 
the  Sherwood  number  and  denoted  by  Sh  in  the  chemical  literature). 

The  object  of  a  calibration  is  to  measure  the  dimensionless  current  i  at  conditions 
where  the  similarity  parameters  on  which  it  depends  have  the  same  values  as  in  the 
environment  where  the  probe  is  to  be  applied;  hopefully,  by  covering  a  sufficiently 
wide  range  of  the  parameters,  the  functional  relationship  for  i  can  be  inferred 
from  the  results. 

In  the  general  case,  the  function  i  obviously  contains  too  many  variables  to  yield 
readily  to  such  an  approach.  Fortunately,  in  many  circumstances  several  of  the 
eight  parameters  in  Eq.  (3)  are  either  of  secondary  importance  or  essentially  in¬ 
variant.  For  example,  Sc  «  1  for  all  gases.  Furthermore,  if  we  confine  our  atten¬ 
tion  to  the  collection  of  ions  at  relatively  high  negative  probe  potentials,  the  electron 
diffusion  coefficient  should  not  influence  the  problem  significantly  (Dj/De  *,0).  The 
parameter  T^/T  was  aiso  assumed  to  be  of  secondary  Importance  (its  value  in  these 
experiments,  where  the  probes  remained  at  essentially  room  temperature,  was  small). 
The  insensitivity  of  i  to  T^/T  was  confirmed  by  the  experiments.  This  is  to  be 
expected  in  conditions  where  the  boundary  layer  (1.  e. ,  the  region  in  which  the  Influence 
of  the  wall  temperature  could  be  important)  is  thinner  than  the  sheath.  This  was  the 
case  in  the  present  experiments,  as  will  be  discussed  below.  Finally,  in  a  low  sub¬ 
sonic  flow  the  ion  collection  would  also  be  insensitive  to  the  Mach  number  and  the 
specific  heat  ratio,  and  the  essential  parameters  would  thus  be  reduced  to  three  in 
number,  1.  e. , 

1  =  i  (Re,  ijj/r,  0)  (4) 

It  is  the  functional  form  of  this  relationship  that  we  sought  to  establish  experimentally. 
It  should  be  pointed  out,  however,  that  although  over  the  major  portion  of  the  wakes 
in  the  ballistics  range  the  flow  was  indeed  subsonic,  the  calibration  experiments 
were  performed  in  a  shock  tube  where  the  Mach  number  M^  seen  by  the  probes 

12 


( 


ranged  from  1. 1  to  1. 4  -  i.  e. ,  It  was  essentially  constant  at  a  low  supersonic  value. 
However,  based  on  existing  analyses  of  supersonic  blunt  probes  (e.  g. ,  References 
6  and  7)  it  appears  reasonable  to  assume  that  at  such  modest  Mach  numbers,  the  flow 
dependence  of  the  ion  collection  is  felt  primarily  through  the  Reynolds  number  rather 
than  the  Mach  number.  The  very  small  range  of  Mach  numbers  available  in  the  shock 
tube  (behind  the  incident  shock)  prevented  a  direct  check  on  this  point.  However,  re¬ 
flected  shock  measurements  (discussed  later)  where  the  Mach  number  was  zero  does 
confirm  this  point.  In  addition,  a  subsequent  application  of  the  calibration  in  the  ballistic 
range  (see  Section  HI),  where  the  Mach  number  ranged  tovery  low  subsonic  values, 
yielded  good  agreement  with  microwave  measurements  of  electron  density  and  thus  the 
assumption  of  Mach  number  Independence  appears  to  be  well  supported. 

The  calibration  procedure  in  the  shock  tube  was  as  follows.  For  given  flow  conditions 
(Re,XD/r)  and  probe  potential  (0),  the  probe  current  (I)  and  electron-ion  number 
density  (n)  were  measured  Independently,  as  described  in  Section  ILA,  and  i  was 
calculated.  The  procedure  was  repeated  for  as  wide  a  range  of  flow  conditions  as 
possible  so  as  to  map  out  the  functional  form  of  Eq.  (4). 

The  pressure,  temperature,  and  flow  velocity  of  the  test  gas  were  calculated  from 
normal  shock-wave  relations.  The  ion  diffusion  coefficient  was  calculated  for  4,000°K 

from  the  cross-section  data  of  Hlrshfelder,  Curtiss,  and  Bird^^  (for  OH.  in  Ar,  it 

2-1  * 
was  9810/p  cm  sec  with  p  the  pressure  in  torr)  and  extrapolated  to  other  tempera- 

(19)  1  75 

tures  with  the  assumed  dependence'  'Dj~T  /p.  A  problem  which  was  encountered 

was  the  determination  of  the  potential  of  the  plasma,  relative  to  which  the  probe  potential 

was  to  be  specified.  The  floating  potential  of  the  probes  in  the  shock  tube  (i.  e. ,  the 

potential  where  the  net  current  was  zero)  was  established  to  be  within  about  1/2  volt 

of  the  shock-tube  ground.  The  plasma  potential  is  somewhat  positive  with  respect  to 
(4-8) 

the  floating  potential,  '  the  difference  being  of  the  order  kT/e  (0. 2  to  0. 5  volt  in  the 
shock-tube  experiments),  the  precise  value  depending  on  the  particular  flow  conditions. 
All  the  voltages  referred  to  in  this  report  were  measured  with  respect  to  the  shock- tube 
ground  on  the  assumption  that  the  difference  between  it  and  the  plasma  potential  was 
small  compared  with  the  applied  voltage. 

Some  results  from  the  calibration  are  shown  in  Figures  6  through  9.  The  approach  was 
to  try  to  infer  from  data  of  this  sort  an  analytic  expression  for  i(Re,X^/  r,0)  which  would 
satisfactorily  correlate  all  the  data. 


Figure  6  shows  typical  data  of  1  as  a  function  of  Re  for  three  values  of  probe  potential, 

-2  -2 

wlthXp/r  In  the  range  4x10  to  9x10  for  all  points.  Data  of  this  sort  suggested  that 
1  may  be  described  by  the  form 

i  =  f(0,  XD/r)  (a+bReft  (5) 


where  f  is  an  increasing  function  of  0,  and  a,  b,  and  0  may  depend  on  Xjy  r.  Figure 
7  shows  measurements  of  i  as  a  function  of  0  at  a  fixed  Reynolds  number  and  Xjy  r. 
This  and  similar  data  suggested  that  the  potential  dependence  could  be  modeled  in  the 
form 


*(X 

f(0,  XD/r)  =  (-0) 


(6) 


where  a  is  a  function  of  X_/r.  (Note  that  0  is  negative.)  a  was  determined  from 
plots  like  the  one  shown  in  Figure  7.  The  results  are  shown  in  Figure  8.  The  points 
are  fitted  with  the  analytic  relation 

a  =  ITOTTSIITFa^  (7) 

which  was  chosen  so  as  to  not  only  fit  the  data  but  to  extrapolate  to  the  expected  theo¬ 
retical  limits  a  -  0  forXp/ r  -«  0  (see  Ref.  20)  and  a  -  1  for  X  ^  r  >  >  1  (see  Refs.  5, 
20)  at  the  lower  Reynolds  numbers. 


Figure  6  Shock  Tube  Calibration,  i  vs  Re  for  Three  Probe  Voltages: 

V  =-2.8  volts  (0=7  to -13),  V  =  8.4  Volts  (0  =  -21  to  38)  and 

V  =  -15  (0  =  -38  to  -70) 
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.SPHERICAL  ION  PROSES  -  DIAMETER  0.7<  mm 


Figure  9  Correlation  of  Shock  Tube  Data  in  Terms  of  Equation  (8) 


It  was  fount*  hat,  to  a  good  approximation,  a  ,  b,  and  fi  in  Eq.  (5)  could  be  re¬ 
garded  as  independent  of  t  n/  r  and  ali  the  results  covering  the  range  of  parameters: 
A  .9  u  2 

4  <  Re  <  10’  ,  10  «;  Xp.  r<  0. 15  ,  -10  <  0  <  -10,  correlated  with  the  relation: 


1 _ 

(1+0.  61 '  n(l  +  r/A_.) 

i  =  (-0)  D 


+  0.  09  Re 


1. 


o" 


(8) 


Figure  9  shows  Eq.  (8)  together  with  all  the  experimental  data  plotted  in  the  same 
form.  Clearly,  the  correlation  appears  to  be  reasonably  successful.  The  data  fits 
with  an  average  scatter  of  less  than  50T  and  shows  no  systematic  deviation  from  the 
correlation  curve. 


Only  incident  shock  measurements  were  used  to  obtain  Eq.  (8)  and  the  value  of  1.  2 
was  obtained  by  extrapolation  of  the  data  to  Re  -  0.  The  reflected  shock  data,  on  the 
other  hand,  is  (in  principle)  capable  of  giving  a  more  accurate  value  since  no  extrapo¬ 
lation  is  required  because  Re  is  zero  (or  very  close  to  it).  Also,  since  M  is  zero, 
its  influence  on  i  can  be  checked. 


For  the  reflected  shock  measurements  the  probe  diameter  was  0.  76  mm  and  bias  -15 
volts.  Since  the  focused  and  electrostatic  probes  were  located  at  slightly  different 
axial  positions  in  the  test  box,  data  obtained  behind  the  reflected  shock  are  the  result 
of  different  gas  histories  behind  the  incident  shock.  Therefore,  there  could  be  no 
point-to-point  comparison  between  the  microwave  and  electrostatic  probes.  Instead, 
the  electron  density  behind  the  reflected  shock  was  determined  as  follows:  The  time 
of  arrival,  At ,  of  the  reflected  shock  behind  the  incident  shock  at  the  electrostatic 
probe  was  determined  from  the  trace  in  Figure  5b.  The  electron  density  behind  the 
reflected  shock  at  the  electrostatic  probe  was  then  calculated  by  multiplying  the 
electron  density  behind  the  incident  shock  at  the  microwave  probe  at  time  At  (see 
Figure  5a)  by  the  density  ratio  ft.  Pg. 

Two  processes  are  operative  around  the  electrostatic  probes  behind  the  reflected  shocks. 

First,  the  flow  velocity  does  not  instantaneously  go  to  zero  behind  the  reflected  shock. 

Fluctuations  caused  by  interactions  of  the  shock  with  the  probe  damp  out  in  a  time  of  the 
2 

order  r  /o  (equal  to  7psec  for  the  conditions  in  Figure  5).  During  this  damping  period 
the  probe  currents  reach  a  minimum  as  Re-0  and  then  increase  due  to  an  increase  in 
ionization  as  a  consequence  of  the  increased  temperature  behind  the  reflected  shock. 

This  increase  is  seen  in  ali  the  traces  (Region  (5)  in  Figures  5a  and  b). 
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Thus  the  method  of  data  reduction  is  not  exact  but  should  yield  a  reasonably  good 
approximation  to  i  for  Re  and  M  equal  to  zero.  Reflected  probe  data  so  treated 
yields  a  value  of  i  =  6. 3  compared  to  5  obtained  from  the  incident  shock  measure¬ 
ments  for  Re  extrapolated  to  zero  (see  Figure  6).  A  corresponding  value  of 
t/flP  "=  2.  2  (with  or  obtained  from  Figure  8  -  i.e. ,  we  assume  to  Re-dependence  of  or) 
was  obtained  behind  the  reflected  shock  compared  to  the  incident  shock's  extrapolated 
value  of  1. 2  (see  Figure  9).  It  is  of  interest  to  note  that  the  former  value  is  closer  to 
the  theoretical  value  (see  next  section)  than  the  latter.  However,  both  are  in  reason¬ 
able  agreement  with  theory.  The  most  valuable  result  of  the  reflected  shock  measure¬ 
ments  is  that  they  confirm  the  Mach  number  dependence  of  i  to  be  small,  as  was 
originally  supposed. 

C.  DISCUSSION  AND  COMPARISON  WITH  AVAILABLE  THEORY 

Since  no  theoretical  results  are  available  for  spherical  probes  for  the  conditions  of 
the  experiments  described  here,  a  direct  comparison  of  Eq.  (8)  with  theory  cannot  be 
made.  It  is  interesting  to  examine,  however,  whether  Eq.  (8)  reduces  to  the  theoretically 
predicted  forms  when  extrapolated  to  the  limits  for  which  explicit  solutions  are  available. 

First,  consider  the  limit  of  a  static  probe,  Re-0.  The  theory  of  Su  and  Lam^2'^  for 
spherical  probes  yields  i  =  2  forAp,  r-0,  and  i  =  (-0)  for  AD  r>  >  1  (for  the  latter 
limit,  see  also  Hoult^  ).  The  correlation  equation,  Eq.  (8),  reduces  to  1  =  1.2  and 
l  =  -1.20,  respectively,  in  these  limits.  This  agreement  is  quite  satisfactory,  particu¬ 
larly  since  the  data  did  not  extend  to  sufficiently  low  Reynolds  numbers  to  warrant 
confidence  in  an  extrapolation  to  Re-0,  and  since  in  any  practical  situation  the  number 
density  in  the  vicinity  of  a  static  probe  (and  hence  the  current  to  it)  will  be  depressed 
from  the  ideal  value  because  of  diffusion  to  the  probe  supports.  The  reflected  shock 
data  gave  1  =  2.  2  and  i  =  -2. 2  0 ,  respectively,  in  the  above  limits.  This  data  does, 
however,  suffer  from  approximations  in  the  data  reduction  process. 

An  extrapolation  of  Eq.  (8)  to  high  Reynolds  numbers,  on  the  other  hand,  does  not  yield 

agreement  with  the  theories  available  for  the  limits  Xn/r-0  and  Xn/r>>  1 .  Lam’s 

(4)  1/2  u  u 

laminar  theory  ,forRe>>l  predicts  i-Re  in  the  limit  Xn/r-0,  whereas  Eq. 

10  ^  (A) 

(8)  extrapolates  to  i  =  0.09  Re  .  For  Xn/r>>  1,  Sonin’s  theory'  '  predicts 

1/2  1/2  u  1/2 
i~  Re  'for  (-0)<<  Re  and  i «  (-0)  for  (-0)>>  Re  '  ,  again  in  disagreement 

with  Eq.  (8). 
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This  disagreement  is  not  surprising.  The  correlation  equation  is  an  approximate  one 
which  was  chosen  to  fit  as  well  as  possible  all  data  obtained  over  the  range  of  condi¬ 
tions  of  these  experiments  (Figure  1).  The  function  i  of  Eq.  (8)  is  actually  a  very 
complex  one  at  high  Re,  as  evidenced  by  the  wealth  of  theories  which  the  electrostatic 
probe  problem  has  spawned,  and  there  is  no  reason  to  expect  that  the  approximate 

« 

correlation  would  not  break  down  outside  the  region  where  it  was  obtained.  For  ex¬ 
ample.  Lam's  theory^  yields  i~  Re1/2  only  in  the  limit  where  the  electrostatic 
sheath  is  thin  compared  with  the  boundary  layer  thickness,  and  Sonin's  theory^  is 
for  the  case  of  very  large  XD<  r  where  the  probe  is  unsheathed.  In  these  experiments, 
on  the  other  hand,  the  sheath  was  large  compared  with  the  boundary  layer  thickness 
except  at  low  Re,  and  comparable,  in  fact,  with  the  probe  radius  (but  never  large 
compared  with  r).  An  approximate  analysis  has  been  given  by  Boer  and  Johnson'  ’  ' 

for  ion  collection  by  a  flat  plate  probe  under  conditions  where  the  sheath  is  thick 
compared  with  the  boundary  layer.  Their  result  can  be  written  in  the  form 

i  -  (XD  r)1/2  <-0)1/2Re3  4.  (9) 

A  quite  similar  result  is  obtained  if  their  approach  is  applied  to  ion  collection  in  the 

stagnation-point  region  of  a  sphere  or  a  cylinder  (see  Appendix  B),  with  the  assumption 

that  the  sheath  is  thick  compared  with  the  boundary  layer  but  thin  compared  with  the  body 

radius  (the  difference  resides  only  in  a  constant,  which  has,  in  any  case,  been  omitted 

from  the  equation  quoted  above).  The  above  equation  does  not  agree  well  with  the  present 

experiments,  where  the  sheath  was  estimated  to  be  not  only  thicker  than  the  boundary 

layer,  but  often  comparable  with  the  probe  radius.  Numerical  computations  for  flat 

(23) 

probes  with  thick  sheaths  have  also  been  performed  by  Dukowicz.  In  a  certain 

range  of  Xn,  r  and  Re,  Dukowicz  correlates  his  computed  results  with  a  formula 
u  1/2  1/2 

which  in  our  notation  reads  i~(-0)  Re  Dukowicz  also  shows  experimental 

0  3  0  7 

data  for  slender  conical  probes  which  he  fits  with  the  relation  i~(*0)  *  Re  *  .  An 
extension  of  his  theory  to  spherical  probes  with  thick  sheaths  would  be  of  interest  in 
the  present  context. 

A  point  of  interest  is  that  in  the  limits  Re>>  1  and  Xp/'r-  0,  our  correlation  equation, 
Eq.  (8),  implies  that: 

I  =  0. 9  en  vf  (v  r2)  (10) 

•  i.  e. ,  that  the  collected  current  is  essentially  equal  to  the  flux  of  ions  "swept  out"  of 

the  flow  by  the  cross  section  of  the  probe.  Noting  next  that  the  flow  velocity  was 
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close  to  the  average  thermal  speed  of  the  ions  at  low  supersonic  Mach  numbers, 
we  see  that  the  same  equation  can  be  rearranged  to  yield 

env  .  9 

I~-y^<4irrZ)  (11) 

One  would,  indeed,  anticipate  a  relation  like  (10)  for  a  probe  in  a  hypersonic,  collision- 
free  flow,  and  a  relation  like  (11)  for  a  probe  in  a  stationary,  collision-free  plasma. 
That  they  both  happen  to  be  implied  simultaneously  by  an  experimental  correlation 

-4 

formula  obtained  at  continuum-flow  conditions  (mean  free  path/probe  radius  from  10 
to  10* *)  at  high  Re  and  M  «  1  should  not,  necessarily,  be  taken  to  mean  that  the 
mechanism  of  ion  collection  in  continuum  flows  is  either  the  "sweeping  out"  of  ions 
from  the  flow  or  their  arrival  at  the  probe  by  simple  random,  collision-free  flux. 

Expressions  of  the  type  shown  in  Equations  (10)  and  (11)  have  been  used'  ’  'to  reduce 

electrostatic  probe  data  in  continuum  flowing  plasmas  whenXp,  r<<  1.  In  hypersonic 

wakes,  however,  the  latter  condition  is  not  met  and,  in  addition,  v{  is  generally  much 

smaller  than  iL.  The  full  calibration  (Equation  (8)  )  then  must  be  used.  The  factor 

(-0)  is,  of  course,  related  to  the  sheath  size,  which  may  be  thought  of  as  increasing 

(26) 

the  collection  area  of  the  probes.  Incidentally,  Renau  and  Scharfman'  'have  pointed 

out  that  our  values  of  (which  range  from  1. 4  to  7)  are  quite  close  (within  25%)  to 
2 

(s/r)  .  Here  s  is  the  sheath  radius  and  may  be  estimated  from  either  the  work  of 
(20)  (27) 

Su  and  Lam  or  Travers  and  Williams.  '  (Both  references  have  expression* 
for  s  which  give  numbers  within  25%  of  each  other. ) 

The  following  section  describes  work  where  the  empirical  calibration  formula,  Equation 
(8),  was  used  to  interpret  spherical  probe  data  taken  in  the  wakes  of  hypersonic  pro¬ 
jectiles.  As  shall  be  seen,  good  agreement  was  obtained  with  microwave  measure¬ 
ments  over  several  orders  of  magnitude  in  electron  density  (despite,  we  note,  the 
difference  in  Mach  number  in  the  two  environments).  This  leaves  us  with  some 
confidence  that  our  calibration  technique  is  successful  and  accurate  to  perhaps  *50%. 
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APPLICATION  TO  HYPERSONIC  WAKES 


An  array  of  twelve  spuerlcal  electrostatic  probes  was  Installed  In  AC-DRL's  ballistics 
range  to  map  electron  and  Ion  density  distributions  behind  hypervelocity  projectiles 
(see  Figure  10).  The  calibration  of  the  previous  section  was  then  used  to  calculate  ion 
density  distributions  in  the  wake  from  the  mean  current  collected  by  the  ion  probes. 

For  purposes  of  computation  we  set  u/D^O.  8  and  rearranged  Equation  (8)  in  the  form 

(12) 

,(28) 


I  =  (-0)®  en  (4 ir  r2)  (1. 2  Dj/r  *  0. 225  vf) 


The  profiles  of  temperature  and  wake  velocity  used  were  those  calculated  by  Wen' 
which  are  in  essential  agreement  with  the  measurements  made  at  CARDe/2^  The 
Debye  length  was  estimated  from  the  electron  density  measured  by  microwave  tech* 


nlques. 


(9,30) 


SHOCK 


Figure  10  Electrostatic  Probes  in  Free -Flight  Range 
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Figure  11  shows  an  example  of  a  radial  ion  density  distribution  obtained  from  probe 

measurements.  Most  of  the  data  at  large  R/D  is  for  5mm  spheres  and  that  at  small 

R/D  is  for  15mm  spheres.  (R  is  the  distance  from  the  probes  to  the  wake  center, 

(31) 

and  D  the  sphere  diameter. )  Since  it  is  observed  that  the  fluid  wake  width'  '  scales 

with  D  and  that  n  D  (the  line  Integral  of  n^  across  the  wake)  is  approximately 
e  p  e 

constant  at  a  given  pD,  the  ion  density  for  the  large  sphere  has  been  multiplied  by 
three  to  allow  a  direct  comparison  in  Figure  11.  The  profile  can,  to  a  good  approxi¬ 
mation,  be  fitted  with  a  Gaussian  curve. 


R/D 

Figure  11  Radial  Distribution  of  Ion  Density  in  Nitrogen  Wake 
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The  ion  density  n(  of  Figure  11  and  for  other  x/D  (ratio  of  distance  behind  model  to 

model  diameter)  has  been  integrated  over  the  wake  diameter  to  give  n^DpCDp  may  be 

thought  of  as  an  equivalent  plasma  "diameter").  These  computations  are  compared  in 

Figure  12  with  the  integrated  electron  density  n  D  as  measured  with  microwave 

(9  27)  e  " 

probes  and  resonators. '  ’  '  Note  that,  although  the  product  n^  changed  only 
modestly  in  the  measurements,  the  corresponding  change  in  n^  was  about  three  orders 
of  magnitude.  The  agreement  is  good  considering  that  the  wake  is  very  turbulent  and 
our  knowledge  of  even  its  mean  properties  (such  as  v^  and  T)  is  very  imperfect. 

The  wakes  studied  in  the  free-flight  range  sometimes  use  mass  addition  to  either  in¬ 
crease  or  decrease  the  electron  concentration.  Under  these  circumstances,  calibrated 
electrostatic  probes  can  be  used  as  an  aid  in  wake  chemistry  studies.  For  example, 
Figure  13  shows  the  positive  ion  density  calculated  from  the  average  probe  current 
measurements  obtained  in  the  wake  of  a  ISmm  sphere  which  injected  freon  into  the 
wake.  Also  shown  is  the  electron  density  nfi  measured  with  microwave  probes. 
Assuming  charge  neutrality,  the  difference  between  these  two  curves  is  the  negative 
ion  concentration  n_  .  These  measurements  indicate  that  attachment  of  electrons 
to  form  negative  ions  was  very  rapid.  On  the  other  hand,  charge  neutralization  between 
ions  was  very  slow. 
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1 — I  I  l  l  l  II 


1 — I  I  I  I  I  III 
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rijDp  (ION  PROBES) 


pD  •  375  (torr.  mm)  Ng 
SPHERE  VELOCITIES  *  6.2  km/sec  H 


Figure  12  Comparison  of  Integrated  Ion  Density  (Electrostatic  Probes) 
and  Integrated  Electron  Density  (Microwave  Probes) 


Figure  13  Application  of  Calibration  to  Wake  Seeded  with  Freon 
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IV.  SUMMARY 


t 


The  present  work  demonstrates  that  the  continuum  electrostatic  probe  is  a  useful 
instrument  for  making  point  measurements  of  charge  density  in  the  wakes  of  hyper¬ 
velocity  projectiles.  Difficulties  encountered  from  lack  of  an  accurate  theory  for 
probe  operation  in  the  conditions  of  a  typical  ballistics  range  were  overcome  by  a 
calibration  of  spherical,  ion-collecting  probes  in  a  shock  tube.  From  this  calibration, 
an  analytic  correlation  formula  (Eq.  (8) )  was  obtained  for  the  probe  current  as  a 
function  of  the  appropriate  similarity  parameters  over  a  wide  range  of  conditions. 

This  empirical  correlation  equation  was  used  to  interpret  the  readings  of  probes 
arranged  in  the  wakes  of  hypersonic  projectiles  in  the  ballistics  range.  Axial  and 
radial  profiles  of  charged  particle  density  were  obtained.  The  results  were  in  good 
agreement  with  microwave  measurements  of  integrated  electron  density. 


% 
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APPENDIX  A 

CALIBRATION  OF  CYLINDRICAL  ELECTROSTATIC  PROBES 
IN  A  CONTINUUM  PLASMA  FLOW 


Most  of  the  work  in  this  report  is  concerned  with  spherical  probes,  which  are  in¬ 
sensitive  to  flow  direction  and  for  that  reason  have  perhaps  the  simplest  geometry 
for  the  measurement  of  ion  density  in  turbulent  flow  conditions.  There  are  also, 
however,  many  applications  for  probes  of  other  geometries  which  are  sensitive  to 
the  flow  direction.  In  this  appendix  are  presented  some  results  from  a  calibration 
of  cylindrical  probes  in  continuum  flow  conditions  simulating  the  typical  existing 
ballistics  range  environment.  Only  ion  collection  is  considered.  The  data  presented 
here  covers  a  fairly  limited  range  of  the  similarity  parameters  which  govern  probe 
response,  and  does  not  represent  a  complete  calibration  of  cylindrical  probes. 
Nevertheless,  it  has  some  fundamental  interest,  and  is  pertinent  to  applications  in 
the  ballistics  range. 

The  experiments  were  performed  with  the  apparatus  described  in  Section  II  in  the 
main  body  of  this  report.  Again,  the  calibration  involved  the  comparison  of  sim¬ 
ultaneous  measurements  of  ion  current  with  probes  and  electron  density  with  a  focused 
microwave  system.  Two  cylindrical  probe  geometries  were  employed.  The  probes 
aligned  with  the  flow  had  a  radius  of  0. 11mm  and  length  of  20mm,  while  those  oriented 
with  their  axes  perpendicular  to  the  flow  had  a  radius  of  0. 15mm  and  length  of  20mm. 

The  probe  calibration  is  again  expressed  in  terms  of  a  dimensionless  probe  current 
i  (which  is  essentially  the  Nusselt  number  for  ion  transfer)(c.  f. ,  Eq.  2) 

1  ■  enfi'Ym  <A-’> 

where  r  is  the  probe  radius,  and  l  is  the  probe  length.  The  dimensionless  current 
is  primarily  a  function  of  four  dimensionless  similarity  parameters,  (c.f. ,  Eq.  4) 

i  =  i(Re,  XD/r,  0,  > /r)  (A-2) 

vf2r 

where  Re  =  — —  is  the  Reynolds  number  based  on  probe  diameter,  Xp/r  is  the 
ratio  of  Debye  length  to  probe  radius,  0  =  eV/kT  is  the  dimensionless  probe  potential, 
and  r/r  is  the  probe  length  to  radius  ratio. 
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Note  that  if  '/r>>  1,  i  will  be  independent  of  * ,  r  for  a  probe  perpendicular  to  the 
flow,  but  not  for  a  probe  aligned  with  the  flow.  This  is  because  the  boundary  layer 
(and  the  sheath)  scales  with  r  for  a  perpendicular  probe,  but  with  '  for  an  aligned 
probe. 

Theoretical  solutions  for  the  dimensionless  probe  current  of  Eq.  (A-2)  are  available, 
or  can  be  derived  by  means  of  available  techniques,  only  for  certain  limiting  cases. 

(4) 

Lam  deals  with  the  case  of  a  thin  laminar  boundary  layer  in  the  limit  XD,  r-0  (sheath 
<<  boundary  layer <<  probe  size),  where  the  ion  collection  is  controlled  by  diffusion 
across  the  boundary  layer  and  is  essentially  independent  of  probe  potential.  For  a 
cylindrical  probe  perpendicular  to  the  flow,  such  an  analysis  yields 

i  ~  Re1/2  (A-3) 

while  for  a  cylinder  aligned  with  the  flow,  one  obtains 

i  -  (f/rf1/2  Re1  2.  (A-4) 

An  approximate  solution  for  the  case  where  the  sheath  is  thick  compared  with  the 

boundary  layer,  but  thin  compared  with  the  probe  dimensions  (boundary  layer  <<  sheath 

(21  221 

<<  probe  size)  has  been  given  by  deBoer  et  al,  ’  '  for  a  flat -plate  probe.  Their 

method  of  solution  is  extended  in  Appendix  B  to  cylindrical  probes  oriented  perpendicular 
to  the  flow  velocity.  An  analytic  solution  is  given  in  Eq.  (B-5)  for  the  dimensionless 
current  i  ,  based  on  the  current  density  at  the  stagnation  point  of  the  cylinder,  as  a 
function  of  Re,  XD/r,  and  0. 

In  the  present  experiments,  the  sheath  was  of  the  order  of  the  probe  radius,  and  hence 
neither  of  the  theories  mentioned  above  are  strictly  applicable.  In  view  of  the  absence 
of  a  more  appropriate  theory,  however,  they  do  offer  a  basis  of  comparison. 

In  our  experiments,  the  Reynolds  number  was  varied  over  three  orders  of  magnitude, 

3 

from  about  1  to  10  .  The  dependence  of  i  onx^/'r  and  0,  however,  was  explored 
only  qualitatively,  since  the  measurements  involved  only  two  values  of  each  of  these 
parameters  for  each  probe  orientation.  The  aspect  ratio  f/r  was  kept  fixed. 

Figures  A- la, b,c,  and  d  show  results  for  the  calibration  of  a  probe  perpendicular  to 
the  flow  direction,  in  the  form  of  a  plot  of  i  versus  Re.  The  data  are  grouped 
according  toX^/r  and  0,  with  Figure  A-la  having  the  larger  values  of  these  quantities 
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Figure  A-l  Calibration  of  Cylindrical  Probe  Perpendicular  to  Flow 
(Radius  r  =  0. 15mm,  Length  i  -  2.0mm) 


and  Figure  A- Id  the  smaller.  The  indicated  spread  in  the  latter  two  quantities  (for 
example  0  =  -40*  15)  does  not  represent  an  inaccuracy  of  measurement,  but  merely 
reflects  the  grouping  of  experimental  data.  Each  set  of  data  points  has  been  fitted 
with  a  simple  relationship  of  the  form  i  ARe*.  The  data  of  Figures  A-la,b,  and  c 
show  that  the  exponent  4  =  1  gives  a  reasonable  fit  to  the  data  just  as  in  the  case  of 
spherical  probes  under  similar  conditions.  The  constant  A  is  0.  6.  0.  5,  and  0.  4  for 
Figures  A-la,b,  and  c  respectively.  The  data  of  Figure  A- Id  is  better  fitted  with  an 
exponent  of  3/4,  however. 

The  data  in  Figure  A-l  is  not  in  good  agreement  with  the  theory  for  sheaths  which 
are  thin  compared  with  the  probe  radius,  as  outlined  in  Appendix  B.  Equation  B-5 
predicts  that  i  (at  the  stagnation  point)  should  be  proportional  to  Re  to  a  power 
between  3/4  and  1/4,  depending  on  0  and  r. 


Figure  A-2  Ratio  of  Sheath  Thickness  to  Probe  Radius  fl/r  for 
Cylindrical  Probe  Perpendicular  to  Flow.  Also 
Boundary  Layer  Thickness 
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The  disagreement  may  be  due  in  large  part  to  the  (act  that  in  the  experiments  the 

sheath  thickness  (calculated  (rom  Eq.  (B-7)  and  shown  in  Figure  A-2)  was  not  small 

compared  with  the  probe  radius,  as  assumed  in  theory.  For  example,  (or  the  case, 

0  =  -40,  and  A_./r  =  0.22,  we  estimate  that  the  sheath  thickness  is  about  one 
u  9 

probe  radius  at  Re~  10  ,  while  (or  0  =  -12.  and  A^/r  =  0. 1  (Figure  A-ld)  it  is  about 

one-hal(  the  probe  radius.  (The  viscous  boundary  layer  thickness  at  the  same  Reynolds 

number  is  about  0.  15  o(  the  probe  radius. )  Note,  however,  that  (or  the  case  with  the 

3  4  1 

thinnest  sheath  (lowest  0  and  A^  r),  i -  Re  rather  than  Re  ,  i.e. ,  the  data  appears 
to  agree  somewhat  better  with  the  thin-sheath  theory  than  in  the  case  of  thicker  sheaths. 
As  expected,  the  current  appears  to  increase  somewhat  with  increasing  0  and  A^.  r, 
although  the  increase  is  small  in  these  cases  and  almost  masked  by  data  scatter. 

Some  results  from  the  calibration  of  probes  aligned  with  the  flow  direction  are  shown 

in  Figure  A-3.  Here,  the  dimensionless  current  is  to  a  good  approximation  propor- 
1  2 

tional  to  Re  ,  as  predicted  by  the  diffusion-controlled  theory  (Eq.  (A-2)  ),  but  there 
is  nevertheless  a  significant  dependence  on  both  A_/r  and  0,  the  current  increasing 
with  both  these  parameters.  The  data  correlates  well  with  the  equation  (see  Figure  A-4). 

i  =0.45  (AD,r)1/3  (-0)1  2  Re1  2  .  (A-5) 

These  results  are  qualitatively  consistent  with  a  situation  where  the  sheath  is  thinner 

1  2 

than  the  boundary- layer  thickness  (leading  to  a  basically  Re  dependence)  but  not 

negligibly  thin  in  comparison  with  it.  For  example,  the  sheath  thickness  has  been 

2 

estimated  to  be  of  about  one  probe  radius,  while  at  Re  =  10  .  the  boundary  layer  thick¬ 
ness  is  a  few  probe  radii  for  an  aligned  probe. 
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Figure  A-3  Calibration  of  Cylindrical  Probe  Aligned  with  Flow 

(Probe  Radius  r  =  0. 11mm,  Length  i  -  20mm,  > ft  =  182.) 
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APPENDIX  B 

AN  APPROXIMATE  THEORY  FOR  ION  COLLECTION  BY  A 
CYLINDRICAL  PROBE  TRANSVERSE  TO  FLOW 


Consider  an  incompressible,  chemically  frozen  plasma  flow  over  a  circular  cylinder 
with  its  axis  perpendicular  to  the  flow  velocity.  The  following  ordering  is  assumed: 
boundary  layer  <<■  sheath  <<  probe  radius.  The  method  of  solution  essentially  follows 
that  of  References  21  and  22.  Referring  to  Figure  B-l,  we  derive  the  following  ex¬ 
pression  for  the  ion  conservation  equation: 

-00 


oo 


iW-en,  A 


•“/K) 


dy 


(B-l) 


where  j  is  the  ion  current  density  at  x,  and  u  is  the  flow  velocity  tangential  to  the 
wall.  For  a  cylinder,  we  use  the  potential  flow  solution  u  =  211^8^10.  The  above 
equation  can  be  rewritten  in  dimensionless  form  as 


i(9)  =  Sc  •  Re 


£  sin  6  •  “6  (9) 


where 


oo 


(B-2) 


is  the  dimensionless  current  density  (Husselt  number),  Re  =  is  the 

Reynolds  number,  Sc  =  u  /  'a  D„  is  the  Schmidt  number  based  on  the  ion  diffusion 

00  00  00 

coefficient  (Scsl  is  assumed  in  the  comparison  with  experiment),  0  =  x/r,  and 
oo 


(B-3) 


is  the  dimensionless  concentration  displacement  thickness.  Following  Reference  3,  we 

construct  a  one-dimensional,  mobility-dominated  model  for  the  sheath  layer  where  n^ 

differs  from  n>  ,  and,  using  Poisson's  equation  and  the  ion  continuity  equation  and 
*oo 

assuming  electrons  to  be  absent  from  the  sheath,  derive  the  following  expression  for  A: 


6  =  cl 


j  +  +3  |0|/ci2)2/3-  |l+3  101/cl2  )1/3 


(B-4) 


where  c  =(XQ/  r)  and  0  is  the  potential  of  the  body  with  respect  to  the  undisturbed  plasma. 
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Figure  B-l  Coordinates  In  Probe  Theory 


Equations  (B-l)  and  (B-4)  can  be  solved  for  the  current  distribution  i(0)  .  For  our 
present  purposes,  it  suffices  to  consider  only  the  expression  for  i  at  the  stagnation 
point,  where  •  =  0: 


3(r/XD)  |0| 


1  4  (r/*D: 


(B-5) 


)  (2/ Sc  Re) 


Note  that  for  (Sc  Re)  *  <<  rAp,  Eq.  (B-5)  simplifies  to 
/  X-  \J/2  /„ _ \3/4 


(B-6) 


Finally,  we  note  for  purposes  of  computation  that  the  actual  sheath  thickness,  i.  e. , 
the  thickness  of  the  layer  past  which  the  probe's  potential  is  shielded  from  the  plasma, 
is  given  in  this  formulation  by  (expressed  for  the  stagnation  point): 


ft  =Ad  (3/4  |0|  )1/2- 


[l  +  (rAD)  (2/Sc  Re) 1/2 
[l  +  (rAD)  (2/Sc Re)1/2 j 


3  1/2 

-  1 


(B-7) 


This  sheath  thickness  is  shown  in  Figure  A-2  as  a  function  of  Re,  0  and  r.  (The 

filled  circles  on  the  left  ordinate  are,  in  descending  order,  the  sheath  thickness  for 

(271 

Re  =  0,  calculated  from  the  expressions  of  Travers  and  Williams'  for  conditions 
(a),  (b),  (c)  and  (d). )  For  the  sake  of  comparison,  the  ratio  of  boundary  layer  thick¬ 
ness  to  probe  radius  at  the  stagnation  point  is  also  shown  in  Figure  A-2.  This  quantity 
-1/2  (32) 

is  given  by  1.  6  Re  ,  if  the  boundary  layer  thickness  ’  is  defined  as  that  for  which 
the  gas  velocity  has  reached  95%  of  the  inviscid  flow  value. 
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IS  ABSTRACT 


An  experimental  shock  tube  program  was  undertaken  in  order  to  calibrate 
electrostatic  probes  which  may  be  used  to  make  point  ion  density  measure¬ 
ments  in  continuum  flowing  plasmas.  An  empirically  derived  formula  which 
relates  the  ion  density  to  the  measured  probe  current  was  obtained  in  terms 
of  the  appropriate  similarity  parameters.  This  empirical  correlation  was 
used  to  interpret  the  ion  currents  collected  by  a  radial  array  of  negatively 
biased  spherical  probes  which  were  positioned  so  as  to  intercept  the  ionized 
wake  of  a  hypervelocity  model  in  a  ballistics  range.  The  integrated  ion 
densities  were  found  to  be  in  good  accord  with  the  integrated  electron  densities 
which  were  simultaneously  measured  using  microwave  methods.  The  applica¬ 
tion  of  the  electrostatic  probe  measuring  technique  to  the  study  of  negative - 
ion  chemistry  in  wakes  is  also  illustrated  and  briefly  discussed. 
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